Higher order correction to the neutrino self-energy in a medium and its astrophysical 

apphcations 



ALBERTO BRAVO GARCIA and SARIRA SAHU 
Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
Circuito Exterior, C. U., A. Postal 70-543, 04510, Mexico D.F., Mexico 

We have calculated the 1/A/* {M the vector boson mass) order correction to the neutrino self- 
energy in a medium. The possible application of this higher order contribution to the neutrino 
effective potential is considered in the context of the Early Universe hot plasma and of the cosmo- 
logical Gamma Ray Burst fireball. We found that, depending on the medium parameters and on 
the neutrino properties (mixing angle and mass square difference) the resonant oscillation of active 
to active neutrinos is possible. 



I. INTRODUCTION 



o 
o 

' The particle propagation in the presence of a finite temperature and density medium has been studied extensively 
, because of its importance in various astrophysical and cosmological scenarios. The vacuum dispersion relation is no 
^ • longer respected in the medium and the properties of the particle are modified. This has been studied in great detail 
[ for photon and charged particles^. The neutrino propagation in a medium was first studied by Wolfcnstcin, who 
' determined the refractive index of Ve from the forward scattering amplitude of i/g-electron scattering^. Using the 
' finite temperature field theory method, the neutrino propagation has been studied in various environments'^"^. 

In the medium, the test neutrino or anti-neutrino will interact with the background particles through both charged 
^ ' current (CC) and neutral current (NC) interactions depending on the nature of the background. For example, the 
, electron neutrino will interact both by CC and NC interaction with the background of normal matter, whereas the 
muon and tau neutrinos will interact only through NC interaction in the normal matter, because of the absence 
of muons and taus in the medium. However, exotic environments like, early Universe hot plasma and core of the 
supcrnovae have a substantial amount of neutrinos, trapped within it, which can contribute to the potential due to 
V ~ V scattering. The leading order contribution to the neutrino potential is proportional to the difference of the 
. number densities of particles. So in an almost CP-symmetric medium (e. g. early Universe) the leading order (1/Af^) 
contribution to the potential is very small and one has to go beyond the leading order (to order 

The early Universe must possess a tiny excess of matter over antimatter i.e. Nb/N^ ~ 5 x 10~^" to avoid total 
annihilation of matter-antimatter and charge neutrality implies such an excess of electrons over positrons, but total 
Oh!. Icpton asymmetry may not be comparable to baryon asymmetry''. In fact larger lepton asymmetry is not being 
^ ■ excluded by any theoretical reasons. As it is well known up to present, neutrino oscillation in the early Universe can 
t~l \ have dramatic consequences on BBN. The active-sterile neutrino oscillation in the early Universe can create a large 
J> ■ asymmetry between Ve and Ve which could have a direct effect on the nuclear reactions and finally the abundances 
^ ' in the light elements®"-'^^. Similar situation may occur, in the cosmological Gamma Ray Burst (GRB) fireball, where 
^ , radiation and electron positron plasma is formed in a compact region with extremely high optical depth, so that, 
' photons can not escape from it and they split into electron positron pairs. This fireball may be contaminated by 
baryons both from the central engine and the surrounding medium^^"^*. Neutrinos from the progenitor as well as 
from the inverse beta decay process within the fireball will propagate in the fireball. Depending on the nature of the 
fireball (e. g. temperature), there may be resonant oscillation of active neutrinos within the fireball^^. 

In this paper we calculate the neutrino effective potential up to order and consider its effect on neutrino 

oscillation in different environments like early Universe hot plasma'^'^ and the GRB fireball. We found that, depending 
on the parameters of the fireball, there can be resonant oscillation from active to active neutrino. On the contrary, 
the anti-neutrino conversion is very much suppressed. 

The paper is organized as follows: In sec. 2 we give a detailed calculation of the neutrino self-energy, from different 
diagrams that contribute to it. Sec. 3 is devoted to the calculation of neutrino effective potential in different back- 
grounds. The oscillation of neutrinos in the early Universe hot plasma and the GRB fireball environment are discussed 
in sec. 4 and a brief conclusion is drawn in Sec. 5. 



II. NEUTRINO SELF ENERGY 



In the medium, the Dirac equation for the neutrino is given by [p— Y,{p)\ip ~ 0, where p is the neutrino four 
momentum and S(p) is its self-energy. The neutrino self-energy can be calculated from the three diagrams given in 
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FIG. 1: One-loop diagrams for the neutrino self-energy in a medium. 



FigUl The contribution due to the W-exchange FigH] (a) is given by 
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where we have defined the four momentum q = k — p, R and L are the right and left projection operators defined as 
R,L= - 2~ ^^"^ {p) is the electron propagator given by 



(p2 — m^) 
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The vector boson propagator up to order is given by 
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In the electron propagator 
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ri[p.u) = e{p.u)f{p.u) + e{-p.u)f{-p.u), 
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where u — (1, 0) is the four velocity of the medium, considered at rest, and f{p.u) is the fermion distribution function. 
The real part of Eq. (|2.ip is given by 
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For simplicity let us define ReYiw — R'^wL and by doing the po integral, we get 
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In general Tjw is a function of the external four momentum = (w,k) and the four velocity u'^ of the heat bath. 
Then this can be expressed as 



— T,w = o-w h + /t, 
where a-w and 6vf &re Lorentz scalars and can be calculated from, 
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-Tr\^ReY,w\ — awi^ + bw- 



For a test neutrino of energy = k.u we obtain, 
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For massless neutrinos fc^ = which implies Ey — |k|. Using Eq. (|2.7p in Eqs. (|2.9p and (|2.10p we obtain 
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with Ee the electron energy and m the electron mass. In Eqs. p.l3p and p.l4p we have three different types of 
terms which are proportional to f(p.u), Eef{p.u) and f(j).u)/Ee and also proportional to corresponding anti-fermion 
distribution functions. The number densities of fcrmions and anti-fcrmions are defined as, 
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respectively, where g = 2 for e,p and n and one for neutrinos. For Ee ^ fi, the fcrmion distribution function can be 
written in terms of a series, as follows 
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where fi is the chemical potential and a; = Pfil. Using Ea. (|2.16p the number density can be written as 



Nf = gm 
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where ai = f3ml and ICi, i — 1,2, .. are the modified Bessel functions of integer order. Similarly, we have 

f d^p f{Ee) gm^ ^ 
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which gives 
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$3 = Js + = 5^ E(-l)''^'^ cosh(az) 
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Inserting these expressions in Eqs. (|2.13p and (|2.14p . wc obtain 
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In a heat bath one has to consider the thermal average of the quantities. So we have to replace E^, by (i?^) 
by (Ee) in the medium. Using the above definitions, Eqs. (|2.13p and (|2.14p can be written as 
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Now for massless neutrino Eqs. (|2.1ip and (|2.12p give, 

Gp 



aw 



-V2. 



A/2 



'E, 



E,{N, - ^e) + 3 {{Ee)N, + {E,)Ne) 



(2.26) 



(2.27) 



5 



and 
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where g'^/AM^ — \plGp, and is the Fermi coupHng constant. Compared to avy is supressed by a factor 
Now we have to calculate the neutrino self-energy due to Z-boson exchange. The Z-coupling to the fermions 
is given by 
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where Cy and denote the vector and axial vector couplings, which are given by 
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The contribution from the Z-exchange diagram Fig[T] (b) , can be obtained from the W-exchange diagram (a) by 
changing Mw Mz- The background here is made of neutrinos and anti- neutrinos instead of electrons and positrons. 
Making the above substitutions we obtain 
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Also az is suppressed by a factor 1/A/| compared to hz- Finally, we have to calculate the tadpole diagram due to 
Z-coupling Fig[l] (c) . The contribution to the self-energy due to this diagram is given by 
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where Zt corresponds to the tadpole diagram due to Z-boson coupling. The real part of this diagram is given by 



(2.34) 
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Finally this gives 
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This result implies that we have to consider the tadpole diagram contribution for different backgrounds obtaining 
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for a background made of electrons, neutrinos, protons and neutrons respectively. For matter containing electrons, 
nucleons and neutrinos, the tadpole contribution will give 
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If the heat bath is charge neutral then, 
and this gives 
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The test particle is moving in a heat bath and the particles (fermions) in the heat bath have average energy 
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where Nf is the fermion number density already defined in Eq. (|2.15p . In the non-relativistic limit (/3m 3> 1), Eqs. (j2.17|) 
to l|2.19p reduces to 
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But for the relativistic fermions ~ |p| and Nf can be written as 



Nf = 



and similarly for photon, this is given by 
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where ^(3) = 1.20206. This gives, N^/N^ = 3/4 and N^/N^ = 3/8. Then proceeding in the same way as in Ea. ^^J^ 
we obtain 



where ^(4) = 7r'*/90 and we obtain the quantities 
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One can define the asymmetry parameter for the particle species a as 
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where a can be lepton or baryon. Let us assume that the average particle and anti-particle energies are the same. 
Then using these we can write 
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Similarly Eqs. (|2.32p and (|2.33p can be written as 
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Finally the tadpole diagram contributions to the matter containing e,p, n and v is 
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For the charge- neutral matter we have L^. = Lp, which gives 
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We have calculated the quantities a and b as given in Eq. (|2.8p for all the three diagrams in Fig[T]and using these one 
can calculate the neutrino effective potential, which is shown in the next section. Also depending on the nature of 
the medium, there can be contributions from different backgrounds to a and b. 



III. NEUTRINO EFFECTIVE POTENTIAL 



We have discussed earlier that, the Dirac equation for the neutrino is modified in the medium and this will give the 
dispersion relation 

det((l + a)/: + 5/t) = 0, (3.1) 

where the a and b are the same as given in Eq. (l2.8|) and det is the determinant. From this equation we obtain, 

(l + a)w + 6 = ±(l + a)|k|. (3.2) 

As we have seen earlier, a is much smaller than 6, so we can neglect the a contribution to the dispersion relation. 
Then the dispersion relation for the particle (positive energy solution) is 

iio-\k\)^^b. (3.3) 

From this we can see that, —6 is the effective potential which will be experienced by the neutrino when propagating 
through a medium. So we have the effective potential 

Veff ^ -b. (3.4) 

For example, let us consider an electron neutrino propagating in a (i) background of only electrons and positrons and 
(ii) only neutrinos and anti-neutrinos. In the first case the contribution to b will be given by 

b = bw + bz,{e), (3.5) 
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and in the second case 



b = bz + hzti^e) + bzAv^i) + bzAi^r)- 
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Let us consider two different physical situations which can be reahzed in the cosmological and astrophysical scenarios, 
the early universe hot plasma, before the BBN^''* and secondly the GRB fireball^'^'^'*'^^. Just before the nucleosynthesis 
when temperature was much above the electron mass but below the muon mass i.e. m <C T <C m^. In this situation, 
the contributions will come from e^, n, p and all the neutrinos. Then the effective potential experience by a electron 
neutrino is given by 
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If the test particle is Vf,, in the above Ea. p.7p . the particle anti-particle asymmetry Lj will be replaced by —Li. 
Secondly, let us consider the GRB fireball as a thermalized radiation and electron-positron plasma of temperature 
about 2 to 10 MeV. We also consider the situation that, the fireball is contaminated by baryons^^^^^'^*. Thus I'e 
propagating through it will experience an effective potential 
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If the fireball is charge neutral i.e. electron asymmetry is the same as the proton asymmetry (L^ = Lp) then 
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The effective potential for and v^- can also be calculated accordingly. 



IV. NEUTRINO OSCILLATION 



Here we consider the oscillation of neutrinos in the above mentioned two scenarios. The evolution equation for 
propagation of neutrinos in the medium is given by 



A cos 26* A sin 261 

A 
2 

where x is active and y can be active or sterile, A = (m| — ra\)/2Ei, = Sm'^ /2E^, V is the potential difference between 
Vx and Vy, (i.e. V = Vx — Vy), E^, is the neutrino energy and 6 is the mixing angle. The conversion probability for 
the above processes is given by, 

A2 sin^ 26 



(F- Acos26i)2-t-A2 sin^2e' 



. ,;(V^- Acos26l)2 + A2sin2 26' 
X sin^ ^ ^ 1 I . (4.2) 



In the above equation, the resonance condition is given by 



(4.3) 
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So the resonance condition crucially depends on the sign of the potential. In the early universe case, if we consider 
the i'^ oscillation, then the effective potential is 



It is obvious from the above equation that the active-active oscillation is independent of the baryonic contribution to 
the potential, whereas, if it is active-sterile oscillation, then the potential does depend on the baryonic component. 
Let us further This assumption is probably reasonable, because, during this epoch all the 

neutrinos are massless and equally populated. If this was the case, then > 6.14 x 10~^Tq (Tg is in units of 10^ eV) 
is needed to have resonant conversion along with the condition Ea. l|4.3p . As already discussed by Enqvist et. al.,^ the 
above oscillation can proceed equally from both sides and then it will not create any asymmetry in electron and muon 
neutrinos. The anti-neutrinos oscillation will be very much suppressed, due to the change in sign of the potential. 
But oscillation of active neutrinos to sterile ones will create asymmetry in v and this oscillation will depend on 
both Icpton asymmetry and the baryon asymmetry in the medium. 

The cosmological gamma ray bursts release about 10^^ erg energy in a few second and make them the most luminous 
object in the universe^'^. The sudden release of a large amount of gamma-ray photons into a compact region, can lead 
to an opaque photon-lepton fireball through the production of electron-positron pairs. It is believed that, the initial 
temperature of the fireball is in the range 3 to 10 MeV, so pair creation takes place and it forms a lepton-photon fireball 
with a typical radius of 10^ — 10* cm. The fireball is contaminated by baryons both from the source as well as from 
the surrounding environment Neutrinos produced within the fireball due to inverse beta decay, or neutrinos 
from the progenitor, can propagate through the fireball. These propagating neutrinos can oscillate resonantly or 
non-resonantly, depending on the nature of the fireball^^. For v^. ^ oscillation, the potential is 



If the asymmetries L^^— L^^ is valid in the early universe plasma then it is exactly the same as in the GRB-fireball. 
So we can say that, GRB can mimic the early universe. One can see that the neutrino potential does not depend on 
the baryon content of the fireball. The reason is, proton and neutron contributions to the effective potential are the 
same for and and because of the difference it cancels out. So even if there are baryons in the fireball, they will 
not affect the v^. <-> or the corresponding anti neutrino oscillation processes. Apart from this we can also see that, 
for the antineutrino process, the potential is negative, hence the resonance condition Eq. (j4.3[) can never be satisfied. 
On the other hand for the Ve v^^ process, depending on the value of Lg, there can be a resonance. First of all, for 
resonance to occur, the condition Lg > 6.14 x lO^^Tg has to satisfy and this is only possible for neutrino processes. 
Secondly one has to explore the neutrino energy, Sm^ and mixing angle of the neutrinos to satisfy the resonance 
condition in Eq. (|4.3p . For example if we take Tg = 3 then for resonance to take place Le > 5.6 x lO^'*. Apart from 
this if we further assume a spherical fireball of radius ~ 100 Km and Lg = Lp then the fireball is contaminated by 
a mass > 1.7 x 1O~^''M0. In the fireball also, (as in the case of Early Universe) the oscillation of active neutrino to 
sterile one depends on the baryon and lepton contain of the fireball. A detail analysis of the neutrino oscillation in 
the GRB fireball is in progress. 



We have calculated the 1 /A/** order contribution to the neutrino self energy in Early Universe hot plasma and the 
GRB fireball, where the higher order contribution to the neutrino effective potential has potential importance. We 
have looked for the condition of resonant conversion/oscillation of active to active neutrinos, which does not depend 
on the baryon content of the propagating medium. On the other hand, the active to sterile oscillation/conversion 
does depend on the baryon content of the medium. It is shown that, under certain circumstances, the GRB fireball 
can behave as Early Universe hot plasma. The resonant oscillation of neutrino in the GRB fireball depends on its 
physical parameters (temperature and particle asymmetry) as well as on the neutrino mixing angle and mass square 
differences. If this conditions are satisfied, the resonant oscillations of active to active neutrinos may have an effect 
on the energy release of the GRB Fireball. This has to be studied in great detail in future. 

This work is partially supported by the grant from PAPIIT-UNAM IN119405. 
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